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excellent agreement with random matrix theory predictions. 
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1. INTRODUCTION 

In this course we shall discuss how random matrix theory (RMT) can be used to describe 
the two main sources of errors in a quantum information process: 

(A) Loss of exactitude due to inherent errors of the physical reproduction of the algo- 
rithm. 

(B) Loss of coherence due to coupling to and entanglement with some outer system. 

We thus distinguish errors in the unitary time evolution from those caused by the loss of 
unitarity due to external action caused by the environment, from which perfect isolation 
is never possible. As we shall see in the third section such considerations are intimately 
related to the possibility of reverting some time evolution by a time reflection of the 
state, i.e. to the old problem of de-equilibration after time reversal, usually associated to 
the name of Loschmidt in his discussions with Boltzmann, about 130 years ago [1]. Yet 
Lord Kelvin [2], some years earlier, gave an excellent account of the problem: 

"If we allowed this equalization to proceed for a certain time, and then re- 
versed the motions of all the molecules, we would observe a disequalization. 
However, if the number of molecules is very large, as it is in a gas, any slight 
deviation from absolute precision in the reversal will greatly shorten the time 
during which disequalization occurs. . . Furthermore, if we take account of the 
fact that no physical system can be completely isolated from its surround- 
ings but is in principle interacting with all other molecules in the universe, 
and if we believe that the number of these latter molecules is infinite, then we 



may conclude that it is impossible for temperature-differences to arise spon- 
taneously. . . " 

Though speaking in a classical context, he clearly points to the two problems affecting 
the stability of a quantum information process: precision and uncontrollable interaction 
with the environment. What is more, we are warned that the separation is artificial, as 
the second implies the first. He further points out the essentially statistical character of 
perturbations and coupling. This indicates why we should use a randomized description. 
Decoherence in terms of spontaneous emission, is the intellectual basis for any model 
based on harmonic oscillators as environment. This assumes a very unusual isolation 
of the system, maybe realistic for the decay of excitations of atoms and molecules in 
interstellar space. 

The use of random matrix theory (RMT) to understand quantum systems started mod- 
estly in the restricted field of low - energy nuclear physics, though it was introduced by 
no one less than Eugene Wigner [3]. Since then, it has evolved to applications reach- 
ing from mesoscopics, molecules, atoms all the way to elementary particles [4]. Early 
proposals of applications to correlation analysis of times series [5] have since developed 
as an important tool from biology [6] to econo-physics [7, 8], and the techniques have 
merged with those developed for quantum systems [4]. 

RMT models for decoherence have been introduced some time ago [9, 10, 11] and 
more recently such a model was proposed for fidelity decay [12] and successfully com- 
pared to experiment [13, 14, 15]. More recently the models and calculation techniques 
have been improved [16] and specialized for the treatment of qubits and quantum regis- 
ters [17, 18]. 

The purpose of this course is to give an introduction to the application of RMT both 
to fidelity decay and to decoherence. We shall start in the next section with a short 
overview of RMT. In section three we shall give an introduction to fidelity decay and 
present the RMT model we intend to use. We then discuss linear response in this context 
and proceed to calculate fidelity for the RMT model in this approximation. A comparison 
with exact solutions and experiments follows and we close with a special case of reduced 
fidelity decay known as the fidelity freeze. In section four we shall concentrate more on 
quantum information systems, choosing as a central system a quantum register of n 
qubits. First we shall discuss some general results that can be derived using the linear 
response approximation, and then perform explicit calculations for a random matrix 
model. To end this section, we give heuristic results for the relation between loss of 
coherence and loss of internal entanglement in the central system. In the concluding 
remarks we dwell on the fact, that chaos can slow down fidelity decay and decoherence. 

2. RANDOM MATRIX THEORY: AN OVERVIEW 

We shall give here a very concise overview of the concepts and definitions we use, and 
some properties of the classical random matrix ensembles as introduced by Cartan [19]. 
For a wider description we refer the reader to the standard textbook by Mehta [20] . A 
modern review is given in [4]. 

To construct an ensemble of random Hamiltonian matrices of minimal information we 



start out with the set of Hermitian matrices and require that the ensemble be independent 
of the basis chosen in the A''-dimensional Hilbert space. This implies, that this measure 
must be invariant under similarity transformations by the unitary group ^{N). The 
measure djU thus fulfills 

d.\x (H) = A^{UHU^) , U e'^iN). (1) 

Checking numbers of conditions, to fix the measure, against the number of variables, 
we find that there are still A'^ missing conditions. These correspond to the eigenvalues 
of the hermitian matrix H, and up to this point are variables, whose distribution is not 
determined. Balian shows [21] that the minimum information requirement is sufficient 
to fix this. His argument leads to matrix elements //„, „ being Gaussian distributed both 
in their real and imaginary part, all with the same width, and independent except for the 
Hermiticity condition, i.e. , 

d/i(H) oc e-«''''(^')/2dRe(//)dIm(//). (2) 

This result can be rewritten in terms of the Haar measure of the symmetry group [in this 
case {N)] and a measure for the eigenenergies Et of H. This idea can be generalized 
to groups other than ^ (N). We thus allow arbitrary invariance groups ^^3. The measure 
then reads as 

d^p{H) oc dhp{G)llexp\E,-E/ UdEi, G G (3) 

¥j I 

where dhp{G) indicates the invariant Haar measure of the group '^^p. 

The alternate groups are the orthogonal group ff{N), labeled with /3 = 1, and the 
unitary symplectic group ^S^S^i^N), labeled with /3 = 4, in which case we get the 
Gaussian orthogonal and the Gaussian symplectic ensembles (GOE, GSE). The former 
is an ensemble of real symmetric Hamiltonians appropriate for describing most time 
reversal invariant systems, except such of semi-integer spin where we have to use the 
latter. If we use the unitary group (/3 = 2) we obtain the Gaussian unitary ensemble 
(GUE), appropriate for describing time-reversal breaking systems. Note that the spectral 
measure is analytic for the GUE and GSE (/3 = 2 and 4) but not for the GUE (j8 = 1); 
this leads to the fact that analytic calculations, in the GOE, are often more complicated 
than for the other two ensembles. 

It is possible to construct the equivalent ensembles of unitary matrices giving rise to 
the circular orthogonal, unitary and symplectic ensembles (COE, CUE and CSE). These 
are unitary matrices with the same invariance properties under similarity transforma- 
tions. Yet they have larger symmetries which define the ensembles completely. Thus the 
CUE is actually invariant under (N) x ^ (A'') as 

dv{S)=dv{USV) (4) 

where U and V run over the group independently. This defines the invariant Haar mea- 
sure of the unitary group, i.e. , the CUE is the unitary group plus its Haar measure. 
The COE on the other hand consists of unitary symmetric matrices that do NOT form a 
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FIGURE 1. The Wigner surmise for the nearest neighbor spacing distribution P{s) and the spectral 
form factor are shown on the left and right panel respectively for the three ensembles discussed. 



group, and its invariance group is the unitary group, but not as a similarity transforma- 
tion. Rather we have 

dv{S)=dv{USU'), (5) 

where ' indicates transposition. While the transposed is not the inverse of a unitary 
transformation, this is the case for the orthogonal subgroup thus displaying the same 
symmetry under similarity transformations as the GOE. The circular ensembles (CUE, 
COE, and CSE) are appropriate to describe unitary evolution operators such as the 
scattering matrix or the Floquet operator corresponding to a periodic time dependence 
of a Hamiltonian. 

Note that for any of the three invariance groups we can construct generalized invariant 
ensembles with arbitrary spectral distributions inserted in eq. (3). Examples include the 
Poisson unitary and orthogonal ensembles (PUE, POE) [22], where a random uniform 
distribution (in an appropriate interval) of eigenvalues is assumed. This ensemble mainly 
serves as contrast to display the effect of spectral correlations. An equidistant "picket 
fence" spectrum can also be inserted to take spectral stiffness to the extreme. 

It thus becomes clear, that the characteristic properties of the classical ensembles are 
found in their spectral statistics, and this has been used extensively. The level density 
is of limited interest, as this is really a specific property of each system. The classical 
ensembles of Hamiltonians lead in, the large A'^ limit, to a semicircular level density, 
which is quite unrealistic. The ensembles of unitary matrices on the other hand lead to 
constant densities of eigenphases on the unit circle, which is often realistic. If analytic 
solutions are thought they refer mostly to the center of the spectrum, where the level 
density is flat and only needs to be normalized to one to study spectral correllations. 
If spectral statistics of some experiment are compared with RMT, the spectrum must 
be unfolded, i.e. the spectrum must be normalized locally. For numerical simulations 
one can take advantage of the generalized invariant ensembles and introduce unfolded 
spectra in the calculations. 

The most popular statistical test is the nearest neighbor spacing distribution i.e. the 
distribution of spacings of neighboring levels in the unfolded spectrum. This function 
depends in principle on all n-point functions, but its short range behavior is dominated 



by two-point function. The Wigner surmise gives the simple estimate 



Puis) 



' '^i^ r2[(/3 + i)/2] 



(6) 



for the nearest neighbor spacing distribution, which is exact for n = 2 [4]. Fig. 1(a) 
shows this surmise which is quite good even for the large A'^. 

Another quantity of considerable interest and that will be used often during this course 
is the form factor 

]_y^t{Ei-Ej) ^ 1 
Nj-; N 



'itE: 



(7) 



i.e. of the two point function in the time domain. Here Ej are the eigenenergies of the 
Hamiltonian under consideration. We define the spectral form factor for an ensemble in 
terms of the ensemble average {Kj) as ^2(^/%) = 1 + 5(?/th) — {K2{t/xY{)) in terms 
of a dimensionless time T = ^/th whose scale is set by the Heisenberg time Th = 
2;r^/(mean level spacing). Exact and approximate analytic results for are known in 
the large A'^ limit. For the GUE the exact form is quite simple and given by 



if < 1, 
if ItI > 1 



(8) 



It is of interest to remark, that for the GSE case, has a singularity at Heisenberg 
time, while for GUE and GOE respectively the first and third derivatives respectively are 
discontinuous. Note finally that, often, K2 is called the spectral form factor. We reserves 
this term for b2- 



3. RANDOM MATRIX THEORY OF FIDELITY DECAY 

As, in this school, the stability of unitary evolution has received less attention, we have 
to introduce the problem in more detail, define fidelity and mention why it is a sensible 
benchmark. We assume an ideal evolution U and a physical evolution Ue. To estimate 
the deviation between these two evolutions at the end of, or at any intermediate time 
(step) of the quantum evolution, a convenient and widely used measure of this deviation 
is the overlap of the relevant register function under the two time evolutions, i.e. their 
cross correlation 

f{t) = {w\M,{t)\w) (9) 
also known as the fidelity amplitude. Here 

Me{t) = Ul{t) Ue{t) = Uo{-t) Ue{t) (10) 

is known as the echo operator, so f{t) = [(i//|i7Q(?)][i7e(?)|v/)]. This name results from 
the fact that we can reinterpret the definition of the fidelity amplitude as an auto- 
correlation function of evolution under the echo-operator, i.e. forward evolution with 
one Hamiltonian and backward (in time) evolution with another one. 



The modulus squared of the fidelity amplitude F{t) = P is the fidelity, a standard 
benchmark in quantum information [23, 24]. Note that, although we do not mark this 
explicitly, / and F are dependent on the choice of the initial state 1 1//) . Often the average 
fidelity over some set of states, typically randomly chosen, is used. The importance of 
this measure resides in the fact that it has minimal bias, as it averages errors over the 
entire Hilbert space. Other choices are rather task-dependent, but certainly important. 
Thus e.g. the success rate for a given task is, in that sense, the most practical. 

After this brief introduction to fidelity decay we shall present the random matrix 
model we use, and then discuss its solution by linear response. Next we shall compare 
exact solutions to the approximate ones we obtained. We shall then mention the concept 
of scattering fidelity to be able to compare theory to experiments from well outside 
quantum information. 



3.1. The random matrix model 

The model describes fidelity decay in a quantum-chaotic (i.e. mixing) system, suffer- 
ing a global static perturbation; that model could be extended to treat also noisy pertur- 
bations. Typically such a perturbation is less harmful, and can furthermore be treated 
using the statistical properties of the noise directly in the correlation functions [16]. 
Chaoticity justifies choosing the unperturbed system from one of the classical ensem- 
bles or more particularly from a GUE. The word "global" implies that, in the eigenbasis 
of the unperturbed Hamiltonian, the perturbation matrix must not be sparse. 

We consider a perturbed Hamiltonian of the form 

//g. = cos(e')//o + sin(e')//i, (11) 

where Hq and Hi are chosen randomly from one of the classical ensembles. Using this 
form rather then the more conventional H^i = Ho + e'H\ has the advantage that the 
perturbation does not change the level density of the Hamiltonian. Changing the level 
density is a rather trivial form of perturbation; even if we only scale the spectrum we 
will have fidelity decay due to dephasing in its most primitive form. When discussing 
the experiments, we shall return to this point. While the use of the classical ensembles 
is fully justified for Hq (by the type of system we discuss), this is not obviously true 
for the perturbation. Non-global perturbations may occur and behave non-generically 
[25] but also global perturbations with different properties could be important, such as 
may result two-body interactions [26] or from maximally time-reversal breaking, i.e. 
hermitian antisymmetric perturbations of a symmetric (GOE) Hamiltonian [27]. We 
shall discuss such situations at the end of this section. 

We are interested in situations where e' scales as 1 / \/N, where denotes the dimen- 
sion of the Hamiltonian matrices. The matrix elements of the perturbation then couple a 
finite number of neighboring eigenstates of the unperturbed system, largely independent 
on A''. Stronger perturbations would practically lead to a loss of fidelity on a time scale of 
the order of the Zenon time given by the inverse of the spectral span of the Hamiltonian. 
Considering large A^, we linearize the trigonometric functions in eq. (11). We further- 
more fix the average level spacing of Hq to be one in the center of the spectrum, and 



require that the off-diagonal matrix elements ofV = H\/ \/N have unit variance. This 
leads us to the conventional form 



He = Ho + eV, 



(12) 



if we rescale the perturbation parameter as £ = y/Ne'. It is easy to check that corrections 
to the Heisenberg time are of order ff{l/N). We may use different ensembles fori^o and 
V . In many cases, the ensemble of perturbations is invariant under the transformations 
that diagonalize Hq. We can then choose Hq to be diagonal and with a spectrum {Ej}. In 
this situation we can unfold the spectrum that defines Hq, to have average level density 
one along the entire spectrum, or we can restrict our considerations to the center of 
the spectrum. This restricts us to situations, where the spectral density may be assumed 
constant over the energy spread of the initial state. Other cases could be important but 
have, to our knowledge, so far not been considered in RMT. 

In this section, the eigenbasis of Hq will be the only preferred basis in contrast to the 
following section, where we deal with entangled subsystems. Unless stated otherwise, 
we consider initial states to be random, but of finite span in the spectrum of Hq. The 
spectral span of the initial state and the spreading width of the //o-eigenstates, in the 
eigenbasis of H^, determine the only additional relevant time scales. They should be 
compared to the Heisenberg time Th, by unfolding the spectrum of Hq. In the limit 
N the Zeno time (of order X]\/N) plays no role. 

The results presented here cover essentially the range from the perturbative up to the 
Fermi golden rule regime [28, 29]. The analysis of the quantum freeze and an exact 
analytical result for the random matrix model will provide additional and/or different 
regimes. The Lyapunov regime [28, 30] as well as the particular behavior of coherent 
states are certainly not within the scope of RMT. 

Finally, we wish to add that, in many situations, the fidelity amplitude / is self 
averaging [eq. (28)]. Therefore we mainly concentrate on the fidelity amplitude, and 
do not bother with the more complicated averages for fidelity F itself . 



We shall follow the approach of [10], which uses the linear response approximation 
expressed in terms of correlation integrals, as introduced by Prosen [31] and discussed in 
[16]. Some preliminaries about units and a brief recapitulation of the interaction picture 
will be useful. 

First we obtain a particularly useful expansion of the echo operator eq. (10), namely 
the Born expansion. The definition of the state ket in the interaction picture is 



where f/e(?) — exp(— i?i/e) is the evolution operator corresponding to Hamiltonian (12). 
Consider the equation of motion of the ket in the interaction picture {h= 1): 



3.2. Linear response theory and RIMT 



W{t))i = t/ • {t)Ue{t) I V/(0)) = Me{t) I V/(0)) 



(13) 




£Vt\\{r{t))i 



(14) 



dt 



where we are using the shorthand 

A, = ul{t)AUo{t) (15) 
for any operator A in the interaction picture. Formal integration of eq. (14) leads to 

|V/(0)/=|v^(0))-ie f Ax%Mx))i. (16) 
Solving the integral by iteration we obtain 

\\\f[t))i=(^-i£ jUx%-£'^ j'^dx j\x%, + ..^\^f{Q)). (17) 
Comparing (13) and (17) we obtain the Bom expansion: 

Me{t) = t-\£ ['dx%-e^ f Ax f Ax%%,+0{e^). (18) 

JO JQ Jo 

From eq. (9) we see that we must only calculate the expectation value of the echo 
operator, and thus of the two objects / (V-c) and // (VtV't')- this point we choose to 
work in the eigenbasis of the unperturbed Hamiltonian where we can write 

Ho = diag(£«), Uo{t) = diag(e-'^«^/'^) . (19) 

We proceed averaging the perturbation over the GUE. This average will be denoted by 
{■)v- As the matrix elements are independent (up to symmetry requirements) complex 
Gaussian variables we have that 

{V,,j)v = 0, {V,,jVkj)v = d,j5j^k- (20) 

This relation fixes a normalization condition on the ensemble. However this normaliza- 
tion condition can always be met, absorbing the appropriate factor in e. 

The expectation value of the linear term of the echo operator yields zero automatically, 
as we can check for the matrix element [■]a.ii'- 

fdxV,) 1 = f'dxo'<^--'^P\v^^p)v = 0. (21) 
The next term in the expansion involves the correlation function (V-rVr')!/- We obtain: 

y 
y 

= 5«^/3i:e'(^r-^«)(^-^'). (22) 



We now average over the initial state. Consider an observable A, and the state = 
with Xi complex random numbers with standard deviation \/N (for normaliza- 
tion), i.e. a random state. Then, denoting average over the initial state as {■)\xj/y 



'J 



'N 



N 



(23) 



In other words, to average expectation values over random states is equivalent to tracing 
the operator. Thus, after evaluating the average over initial conditions of the correlation 
function, we obtain 



v,\w) N 



-trVe^' 



1 



£e'(^r-£«)(t-T')=^2(T-T'). (24) 



a, 7 



We have thus related fidelity decay directly to the form factor of the environment. For 
the average fidelity amplitude we obtain 



(/(r)) = l-£2 TdT l\x'K2{x-x') 
Jo Jo 



2 2 



t ft 

dT / dT'Z72 

Jo 



T-T' 



(25) 



Any stationary ensemble from which Hq may be chosen, yields a particular two-point 
function b2- The correlation integral over b2 for the GOE is discussed in [10]. If we use 
a PUE, i.e. a random level sequence, ^2(0 = 0^ the last term in (25) vanishes. Typically 
(at least in the case of the classical ensembles), spectral correlations lead to a positive 
b2, such that fidelity decay will be slowed down by these correlations. However the 
dominant term is, before/after the Heisenberg time, the linear/quadratic one respectively. 
The correlation integral for the GUE is given by 



dr / Ax'b2 



min{?, Th} 



min{r,TH}' 
3th 



(26) 



Averaging the perturbation over other ensembles is also possible. In the final result 
eq. (25), on has to substitute ?^/2 with t"^ / jSy where labels the ensemble from which 

V is drawn [10]. The result (25) shows two remarkable features: The first is that the 
linear and the quadratic term in t both scale with The second is about the role of 
the two possibly different ensembles used for the perturbation and for Hq. As long as 
the invariance group of the perturbation allows to diagonalize Hq the characteristics of 

V affect only the prefactor of the ?^-term, while the characteristics of Hq affect only the 
two-point form factor Z?2. 

The expansion in time, eq. (25), contains the leading terms for both the regime known 
as perturbative and the one known as the Fermi golden rule regime [28]. Both are 
exponentials of the corresponding terms in the linear response approximation. It is then 
tempting to simply exponentiate the entire £^-term to obtain 



(/(0)=exp<^-e 



2 +/3v 



dT / Ax'b2 
Jo 



(27) 



This expression will prove to be extremely accurate for perturbation strengths up to the 
Fermi golden rule regime. Some justification for the exponentiation is given in [32]. 
While exponentiation in the perturbative regime is trivially justified, our result shows 
that for times t <^ Th, we always need the linear term in t to obtain the correct answer; 
with other words, the separation into the two regimes is somewhat artificial: Before the 
Heisenberg time, the linear term is dominant, and except if fidelity is way down before 
Heisenberg time, after this time the quadratic term will dominate. Thus the RMT formula 
is adequate to describe the entire transition between the perturbative and the Fermi 
golden rule regime. In experiments the interplay of both terms is important [15, 33, 14] 
as we shall see below. On the other hand, for stronger perturbations fidelity has decayed 
before Heisenberg time to levels where our approximation fails. Comparison with the 
exact result [34, 35] will show that the exponentiation allows to extend the linear 
response result from a validity of (/(?)) ~ 1 to a validity range of (/(?)) ^0.1 

Note that the pure linear response result is probably all we need for quantum infor- 
mation purposes, as processes with fidelity less than 1—7], where rj ~ 10^^, are not 
amenable to quantum error correction schemes [23]. The exact treatment will show 
where to expect additional effects, but experiments at this time are still limited to 
(/(0)>0.1[33, 14,15]. 

We now turn our attention to fidelity. It can be calculated in the linear response 
approximation along the same lines as above. One obtains [10]: 

{F{t)) = (1/(0 12) = imf + e' (2//3k) iprt^ +ff{e') . (28) 

Here ipr = Lv K^vlv^)!"^ indicates the inverse participation ratio of the initial state 
expanded in the eigenbasis of Hq. We deduce two extreme effects: On the one hand it 
shows the self-averaging properties of this system. For states with a large spectral span 
in Hq, the correction term that marks the difference between {F{t)) and | (/(?)) P goes to 
zero as the inverse participation ratio becomes small (~ l/N). On the other hand, for an 
eigenstate of Hq, ipr = 1, and hence the quadratic term in eq. (28) disappears. Moreover, 
the correlations cancel the linear term after the Heisenberg time. Thus, we find that 
after Heisenberg time the decay stops for an Hq taken from a GUE and continues only 
logarithmically for a GOE [10] up to the onset of the next term in the Born expansion. 



3.3. Supersymmetric results for the fidelity amplitude 

The exponentiated linear response formula (27) agrees very well with dynamical mod- 
els [10, 36] and experiments [33, 14, 15]. Yet this result has to fail for sufficiently strong 
perturbations, even if we forget its heuristic justification. Fortunately many problems in 
RMT can been solved exactly, and recently Stockmann and Schafer [34, 35] have done 
exactly this for the model given in eq. (12), for GOE or GUE matrices, in the limit of 
infinite dimensions. The solution for the GSE will be published shortly [37]. 

More specifically, they choose Hq and V independently but both from the same 
classical ensemble, and compute the fidelity amplitude (/(?)) with the help of super- 
symmetry techniques. For a detailed account we refer to the original paper [34]. They 
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FIGURE 2. Fidelity amplitude decay for the random matrix model, defined in eq. (12) (taken 
from [34]). Part (a) shows (/(f)) for the GOE case, as obtained from the exact expression, eq. (30), 
(soUd lines), together with the same quantity, as obtained from the exponentiated linear response result, 
eq. (27) (dashed lines). The perturbation strength has been set to the following values: ~ 0.2, 1,2,4 
and 10. Part (b) shows (/(f)) with = 100 for the GUE case (solid line) and the GOE case (dashed line), 
as obtained from the exact expressions, eqs. (29) and (30), respectively. In both figures the authors use the 
convention T = f and /^(t) = (/(f))- 



obtain 

1 /■min(f,l) , 

{f{t)) = - du{l+t-2u)c-'^^^+'-^"^'/^ (29) 

t Jo 

for the GUE case and 



(/(0) = 2/ Aufdv 

Jmax(0,f-1) Jo 



/max(0,f-l) Jo [t'^-V^Y ^J {u^ - v'^){{u + \Y - v'^) 



for the GOE case (Th = 1 in this section). These solutions are valid for arbitrary time 
independent perturbation strength, in the limit N ^ oo. 

In Fig. 2, we reproduce two graphs from [34]. The left one, compares the exact 
and the exponentiated linear response result for {f{t)) for the GOE case. For large 
perturbations we find a qualitative difference in the shape of fidelity decay as a shoulder 
is forming in the exact results. For even stronger perturbations, depicted on the right 
panel, this becomes more pronounced as a revival is seen at Heisenberg time. Yet, the 
revival is noticeable only for very small fidelities, of the order of 10~^ for the GUE 
and 10^^ for the GOE. In all cases agreement with the exponentiated linear response 
formula is limited to (/(?)) ^0.1. Exponentiated linear response is thus adequate for 
most applications. Indeed it was difficult to come up with a dynamical model which can 
show the revival. In [38] a kicked Ising spin chain [39] has been used to illustrate the 
partial revival, as shown in Fig. 3. Taking advantage of the relative ease of numerical 
calculations in such a model, the authors used a multiply kicked Ising spin chain in a 
Hilbert space spanned by 20 qubits, and averaged over a few initial conditions. The aim 
was to obtain the partial revival with as little averaging as possible, relying on the self- 
averaging properties of the fidelity amplitude. The model shows random matrix behavior 
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FIGURE 3. (Color online) The decay of the fidelity amplitude in a dynamical system, in the case of 
broken time-reversal symmetry, in a regime where the partial revival is observable (taken from Ref. [38]). 
More precisely, logjo(/(f )) plotted. The triangles show the real part of the fidelity amplitude as obtained 
from the numerical simulation; the circles show the imaginary part, which goes to zero as 1 / \/N due to 
state averaging. The perturbation strength is = 3 1 .78. The solid line shows the exact theoretical result, 
eq. (29), the dashed line shows the exponentiated linear response result, eq. (27). 



up to small deviations as far as its spectral statistics is concerned. The result for the decay 
of the fidelity amplitude is reproduced in Fig. 3. Yet it will probably be difficult to see 
this effect in an experiment. 

Thus the random matrix model captures all fidelity decay scenarios from Gaussian 
(or perturbative) to exponential (or Fermi-Golden rule). Furthermore it displays an 
additional feature, namely the weak revival at Heisenberg time. Naturally the random 
matrix model does not capture the Lyapunov regime. This is an additional semi classical 
regime for fairly strong perturbations [30, 16] where the decay is determined by the 
classical Lyapunov exponent and in some range becomes independent of the perturbation 
strength. It is not obvious that such a regime always exists. For example, for the kicked 
spin chain this regime has not been observed, and it may well be that there is no semi- 
classical limit for this system. 

The revival at Heisenberg time can be related to a revival in the cross -correlation 
function of spectra in the time-domain through a recently discovered direct relation 
between this cross-correlation and fidelity decay [37]. 

3.4. A micro wave experiment of fidelity decay 

The dynamics of classical electromagnetic waves in a thin resonator is equivalent to 
the Schrodinger equation of a single quantum particle with two degrees of freedom. This 
is exploited in the microwave cavity experiments, where properties of two dimensional 
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FIGURE 4. Geometry of the billiards (figure taken from [33]). In the billiard in the right figure bouncing 
ball orbits have been avoided by inserting additional elements. 



quantum billiards can be studied. The Marburg group considered cavities (Fig. 4) con- 
sisting of a rectangle with inserts that assure chaotic, but not necessarily hyperbolic, 
ray behavior. Both situations, with and without parabolic manifolds leading to so-called 
bouncing ball states, were considered. The presence of bouncing ball states is known to 
lead to a weak deviation from RMT behavior [40]. This will serve to check how impor- 
tant small deviation, from the "universal" RMT behaviour, can be for fidelity decay. One 
of the walls was movable in small steps allowing the perturbation Y occurring in an echo 
experiment [see eq. (12)]. It is important to note that the shift of the wall changes the 
mean level density and by consequence the Heisenberg time. This trivial perturbation, 
which would cause very rapid fidelity decay, is eliminated in this case by measuring all 
times in proper dimensionless units, i.e. in terms of the Heisenberg time. Two antennae 
allow to measure both reflection and transmission channels. The experiment was car- 
ried out in the frequency domain. Afterwards the Fourier transform is used to obtain 
correlation functions in the time domain. We thus get autocorrelation functions for any 
given configuration, as well as cross correlation functions associated with the different 
positions of the moving wall, i.e. between evolutions with different Hamiltonians. After 
normalizing the latter with the autocorrelation we find the so called scattering fidelity. 
For chaotic systems, in the weak perturbation limit, scattering fidelity has been shown 
to be equivalent to the fidelity amplitude [33]. Figure 5 shows the experimental and 
theoretical result in excellent agreement in the absence of bouncing ball states and, as 
expected, in lesser agreement if they are present. In principle in our model the pertur- 
bation strength £ is a free parameter. Yet in this case it was obtained independently by 
exploring the level-dynamics, i.e. the movement of the energy levels under change of 
the Hamiltonian. The same random matrix model can account for this situation. For low 
frequency ranges where resonances are well separated, this so-called level dynamics can 
be measured. The perturbation strength was extracted and used in the formula for fidelity 
decay. 

Similar experiments have been performed by Weaver and Lobkis [41] on the coda of 
vibrations of metal blocks. Using the concept of scattering fidelity, these experiments can 
be reinterpreted as fidelity measurements [15]. Quantum optics techniques in principle 
allow direct measurement of the fidelity amplitude, using a single qubit as a probe. The 
general context in which this is possible is shown in [42] and a detailed proposal for 
the experiment is given in [36] for the case of a kicked rotor where the kicking strength 
is perturbed. Actual experiments where not performed yet, though an experiment with 




FIGURE 5. Logarithmic plot of the scattering fidelity fii{t), directly related to the fidelity amplitudes 
(figure taken from [33]). The smooth curve shows the linear-response result. For the billiard without 
bouncing balls, the perturbation parameter e was obtained from the variance of the level velocities; in the 
other case it was fitted to the experimental curve. 



four different Hamiltonians, two for the forward and two for the backward evolution was 
performed [43]. 



3.5. Residual perturbations and the Quantum freeze 

After the Heisenberg time, fidelity decay is essentially Gaussian, if it has not yet 
decayed significantly [see eq. (27)]. This decay is determined by the diagonal elements 
of the perturbation (in the interaction picture) in the eigenbasis of Hq. If this term is zero 
or very small, we should see a considerable slowing down of fidelity decay. That such a 
possibility exists was first noted for particular integrable [44] and chaotic [45] dynamical 
systems. Yet the generality of the phenomenon is best understood in the context of the 
RMT model discussed in [16, 27]. 

Sticking to the representation where Hq is diagonal, we shall thus consider situations 
where the perturbation is some random matrix with zero diagonal. The off diagonal 
elements will form the so called residual interaction Vres = V- Many ways to implement 
such a residual perturbation are conceivable. Most are basis dependent. The case, where 
Hq is chosen from a GOE and V is an antisymmetric Hermitian random Gaussian matrix 
is of interest, because the real matrix Hq can be diagonalized without disturbing the 
residual (i.e. zero diagonal) character of the perturbation. Thus it has strong invariance 
properties and indeed this case is the only random matrix model, so far, for which an 
exact analytical result has been obtained [27, 46] by supersymmetric techniques. 

Consider the case of a random perturbation (of GUE or GOE type) with diagonal 
elements set to zero. We have a clearcut stabilization of fidelity after the Heisenberg 
time because the double integral in the linear response expression cancels the linear 
term in t exactly for the GUE case and (up to a logarithmic correction) the GOE case. As 
mentioned, above the quadratic term is absent, and decay, as obtained by linear response, 
essentially ends at the Heisenberg time. Decay at much larger times results from the next 
term in the Born series and will go as 
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FIGURE 6. The complement of the average fidelity amplitude 1 — {f{t)), for a weak perturbation, 
e = 5 X 10^^. Frame (a) shows the GOE case with a purely imaginary antisymmetric random Gaus- 
sian perturbation. Frame (b) shows the GUE case, with an independent GUE with deleted diagonal as 
perturbation. The exponentiated linear response approximations are plotted with short dashed lines (for 
comparison, the exact results for full GOE and GUE perturbations are shown with chain curves). The 
forth order results [33] from time-independent perturbation theory are plotted with long dashed lines. In 
the GOE case, the exact analytical result is plotted with a solid line. The random matrix simulations are 
plotted with points. (Taken from [27].) 



We see in Figure 6 that IVIonte Carlo calculations agree perfectly with the exact 
solution for the GOE case with antisymmetric perturbation. Furthermore both the GOE 
and the GUE cases are well described by the linear response result up to the plateau and 
the end of the plateau agrees with the behavior of £^ term. Experiments are not available 
yet, but good agreement was shown with a dynamical model [16, 27]. 

The case with zeros forced on the diagonal of the perturbation is not as artificial 
as it might seem, because is is a standard practice of mean field theory to include the 
diagonal part of the residual interaction, which does not affect eigenfunctions, in the zero 
order mean field Hamiltonian. This leaves the residual interaction with zero diagonal. 
Unfortunately, in such a situation, even if the single particle spectrum is of RMT type, the 
spectrum of the n-particle problem without residual interaction looks more like a random 
spectrum, and this gets worth if the diagonal term of the interaction is included in Hq. 
In such a case the linear term after Heisenberg time is not canceled, and we find linear 
decay after Heisenberg time [16]. In [26] fidelity decay of two-body random ensembles 
for n Fermions were studied including the diagonal part of the two-body interaction in 



Hq. The linear decay for the average fidelity decay was found. Yet remarkably it was 
shown, that the median behavior does display the freeze. This means, that in a typical 
realization the freeze is present and thus the time evolution of a mean field theory with 
weak residual interaction will follow the exact one, after some initial deviation, for a 
long time. 

4. PURITY 

We now turn our attention to the second problem addressed in the introduction. We 
shall study the loss of coherence due to coupling to, and entanglement with, some outer 
system. This is what people often refer to as decoherence. Quantifying and eventually 
controlling the degree of decoherence is a mayor problem of quantum information. 
Assuming a pure state in a bipartite system = M'a ® M'b, one can quantify the degree 
of entanglement via the purity. Let 1 1//) G M' be the pure state in the whole Hilbert space. 
The reduced density matrix of system A (B) is then = tr^ | V^)(v| (Pb = tr^ |v/)(v/|). 
Purity is simply 

P = trp2 = trp| (31) 

which measures the degree of mixedness of each reduced density matrix. One could also 
use other mixedness measures such as the von Neumann entropy. We prefer purity above 
others because it has an analytic structure that simplifies analytic treatment. In particular 
it is not necessary to diagonalize the density matrix to obtain purity. 

In this section we discuss purity decay, thus quantifying decoherence. First we derive 
a general expression for purity decay of n non interacting qubits. We shall obtain a sum 
rule where each term involves a single qubit and its entanglement with the remaining 
ones [47]. After that, we assume random matrix environments and couplings, following 
the scheme proposed in [10, 17] and developed in detail in [18]. We shall not develop 
the model with full generality to keep the discussion at a basic level, allowing us to go 
slowly through the simplest (though representative) cases. 

4.1. n qubit purity decay 

First we derive a general formula for purity decay of n qubits under very general 
assumptions. The problem is solved using the enlarged Hilbert space 

jr = J^(g)J^ (32) 

where is the Hilbert space of the central system and ^ (with dimension A'e) 
is the Hilbert space of the environment. The central system itself has a structure as 
it is composed of several qubits: = <S>'i=i^i with dimJ^ = 2. The Hamiltonian 
governing the system has the usual structure 

H^=H,+H, + XV, (33) 

where He acts on the central system, on the environment and is the coupling be- 
tween central system and environment. We shall moreover require that the Hamiltonian 



of the central system is local, i.e. He = LLi Hi where Hi acts on J^. We further require 
that each qubit is separately coupled to the environment: 



(34) 



with V^'^ acting on ^ (g) We do not restrict any component of the Hamiltonian to 
be time independent, though we shall not explicitly show the time orderings required if 
such dependence exists. 

The state in the whole Hilbert space after a time t is \ ^^{t)) = Ux{t)\\ifQ). Uiit) is the 
evolution operator at time t associated with the Hamiltonian (33) and | i/Zo) = | V^c) ® \ Ve) 
(IVc) e '^c, IV^e) e ^) is a separable pure initial state. We write this state as |i//o) = 
X^^jc^l/i) or, using tensor notation, |i//o) = T.ij.k^ijk\ijk) with \ijk) being an element of 
an orthonormal separable basis in = (g) Jifq ® M's- We shall keep the convention 
that Greek indices are used in the whole Hilbert space whereas Latin ones are used for 
subsystems. 

Purity is the trace of the square of the density matrix of the system in question 
[eq. (31)]. Thus, we are interested in calculating 



P{t) = tr, (tre I Wit) ){¥{f)\f= trc tr.M^ (t) \ Vo) iVo I^J (t) 



(35) 



where (t) is the echo operator defined in conjunction with fidelity, see section 3. The 
last equality is obtained noticing that the difference between the forward evolution oper- 
ator Ux{t) and the echo operator M;(^(?) is the local (with respect to environment-central 
system separation) operation i7j(?)- Since entanglement properties remain unchanged 
under local operations, purity has the same value for the state evolved with either the 
forward or the echo evolution operator. Since for large purities Mp^{t) ~ 1, a series ex- 
pansion is again feasible for long times. We obtain 



(36) 



(37a) 
(37b) 
(37c) 



P{t) = l-2A2^'dT^'dT'ReA(T,T') + 0(A^), 
where A(t, t') = Ay(T, t') -Ai(t, t') +A2(t, t') -^3(1, t') and 

Ay(T,TO = p[%%'pO®po\ = X^X*>j,^iffka*j,,^[Vr%']ijk,^i, 

A2 ( T, t') = 77 [Vt'PoK ® Po] = X^X*yXi> fk'X*fk' [^t] iik,il [Vt'] }jk,v ' 

^3 ( t') = j!? [Vt'Po ® PoVt] = X^X*,jj^i, fk,X*y [Vr] ijk,^ [V^f] *j,i^, y 

(summation over repeated indices is assumed in these equations). The bilinear functional 
p[pi ®p2] = tr(trePitreP2) has been introduced, together with po = |v/b)(V^o|, to sim- 
plify the expressions. Notice that each term depends on the coupling in the interaction 
picture, see eq. (15). 

Taking advantage of the particular structure of the coupling [eq. (34)] we can rewrite 
the integrand in eq. (36) as a sum of terms, each with two indices labeling the qubits. 

X^A{t,t') = £ XiXjA^''j\T,T') 



with 



+ p[vi'^po ® V^Ppo] - p[V^>Po ® PoVi'^] . (38) 

The terms a(''')(t, r') are cross correlation functions for i ^ j and autocorrelation 
functions for i = j. If the cross correlation functions drop to zero fast enough, these 
terms can be neglected. This happens if the environment Hamiltonian is chaotic [47] or if 
the couplings are independent from the outset. Assuming fast decay of cross correlation 
functions [A^'^-''-') ^ 0] we obtain 

^W = l-I(l-4^W), (39) 

;=1 ^ ^ 

Pi;^ (0 = 1- dzj^ dT'A(''') (t, t') + O (A^) . 

Each term Psp (t) represents purity decay of the whole register when only a single 
qubit is coupled to the environment. That configuration is called spectator configuration. 
Eq. (39) is a sum rule for decoherence and is the central result of [47]. 



4.2. RIVIT in the spectator configuration 

We now solve the problem of a single qubit in the spectator configuration, when the 
Hamiltonian of the environment and the coupling are chosen from the classical ensem- 
bles. We use the following subscripts: "c" (for central system), "e" (for environment) , 
"q" (for the coupled qubit), and "s" (for the spectator). 

We have seen that the problem can be formulated in terms of the echo operator. 
However, in the spectator configuration, the echo operator does not contain the internal 
Hamiltonian of the spectator space. Effectively one can thus drop those terms and write 
Hq = Hg+Hq (with Hq the internal Hamiltonian of the coupled qubit) and Hx =Hq + XV 
with XV the coupling between the coupled qubit and the environment. Notice how no 
part of the Hamiltonian is acting on the spectator space. 

At this point we rederive a result from [18]. Purity decay in the spectator configuration 
depends only on Pq = trg | i//c) ( V^c| the reduced density matrix of the coupled qubit alone: 



Pit) 



trc 
tre 
tre 
tre 



tr^Mx0is\Vo){Vo\Ml0l, 
trcM;L®ls|V^o)(V^o|^l®as 



trqMA(trs|v/o)(V^o|)Mj 



trqM;L(Pq® ke)(V^e|)M' 



(40) 
(41) 
(42) 
(43) 



The first step uses only that the dynamics in the spectator space is trivial. In the second 
step we use the property that for a pure bipartite system, purity does not depend on 
the subsystem observed. Thus instead of measuring purity of the quantum register we 
formally calculate purity of the environment. In the next step we split the trace over all 
the register into the trace over the coupled qubit and the trace over the spectator qubits. 
Next, we explicitly shift the trace over the spectator space to one over the initial state 
obtaining the final result: purity decay in the spectator configuration will only depend 
on the reduced density matrix of the coupled qubit. 

Since to purify a qubit we only need one additional qubit, we can replace ^ by the 
Hilbert space of a single qubit thus greatly reducing the complexity of the problem! We 
can now refer to the solution published in [18] for two qubits without loosing generality. 
We carry on the explicit calculations, in a simple case. 

We assume that the nontrivial part of V is well described by the GUE, with a normal- 
ization set by the condition 

{V,j.ki) = 0, {V,j,kiV,>f,k'i') = (5;./5/,/0 {5i^k'5i',k) (44) 

where the is and ks run over the environment and the coupled qubit whereas the is and js 
run over the spectator qubit. The first set of 5s describe the fact that the spectator is not 
coupled, and the second set of 5s is the usual GUE relation. The normalization condition 
is no restriction as an arbitrary factor can be absorbed in the coupling constant A. 
However this normalization guarantees that decoherence decay is roughly independent 
on the environment size, for big environments. One could also assume that V is well 
described by the GOE, however this leads to complications (due to the weaker invariance 
properties of the ensemble) on which we do not want to dwell [48]. We first average over 
the coupling. To do so we chose a basis that diagonalizes Hq. In this basis, [V-r]//,/'/ = 
^jj'^^'^'^ij'^^^'^'' ■ The is and 7s run over ^ ® ^ and Jifs respectively, and both Ei and 
£■,' are eigenenergies of Hq. As long as V and Hq are independent, relations (44) remain 
unchanged, so we can write 

2A'e 2 

(45) 

and 

This greatly simplifies eqs. (37). 

A further step towards obtaining the final expression is to average over the state 
of the environment. In order to ease the following work, we rewrite the separability 
condition as Xjjk = Wi<^jk where |v/c) = L/ V^/IO the initial state of the environment 
and |i//c) = Y.j.k^ik\jk) the state of the two qubits. However, we wish to write |v/c) in 
its simplest form. We use the invariance properties of the ensemble defined by Hi, but 
under the restriction that Hq is still diagonal [to be able to write eqs. (45) and (46)]. This 




FIGURE 7. (Color online) We present a plot to help visualize the way the initial condition is 
parametrized. On the left, the coupled qubit. The eigenvectors of its internal Hamiltonian (|0) and |1)) 
are represented in blue. The z axis is chosen parallel to the vector |0). The x axis is chosen in such a way 
as to make both |0i) and have real coefficients i.e. such that the xz plane contains |0i) and On 
the right we represent the second qubit where we have absolute freedom to choose the basis (even if an 
internal Hamiltonian is present), and thus we choose it according to the natural Schmidt decomposition. 



freedom allows any unitary operation on the purifying qubit, and a phase transformation 
exp(ia//q) of the eigenvectors of //q (a G M). Taking into account this freedom we now 
show how to construct a basis { |0) , 1 1 ) } ® { |0) , 1 1 ) } in which an arbitrary initial state of 
the central system can be written as 

IV/c) =cos0i(cos02|O)+sin02|l))|O)+sin0i(sin02|O) -cos02|l))|l), (47) 

and still, //q = ^ |0) (0| — 1 1 1 ) ( 1 1 is diagonal. To find this basis we start using the Schmidt 
decomposition to write 

I v/c) = cos di IO1O2) + sin 01 1 1 1 12) (48) 

with {|0,), |1,)} being an orthonormal basis of particle For the coupled qubit, we fix 
the z axis of the Bloch sphere (containing both |0) and 1 1)) parallel to the eigenvectors 
of //q, and the y axis perpendicular (in the Bloch sphere) to both the z axis and |0i). 
The states contained in the xz plane are then real superpositions of |0) and |1), which 
implies that |0i) = cos02|O) + sin02|l) and |li) = sin02|O) — cos02|l) for some 02- 
In the second qubit it is enough to set |0) = IO2) and |1) = |i2). This freedom is also 
related to the fact that purity only depends on Pq = trs | i//c) (V^d- A visualization of this 
procedure is found in fig. 7. The angle 0i G [0, ;r/4] measures the entanglement between 
the coupled qubit and the spectator space whereas the angle 62 G [0, Ti/2] is related to an 
initial magnetization of the coupled qubit. 



We have then that 
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2 Reis:2(T-T0 (49) 
-ReiS:2(T-T0 [1+cos(A(t-t')] (50) 

Several approximations have been done. In the first step we assume that the imaginary 
contribution from the environment is very small. This is justified noticing that the form 
factor of a random matrix ensemble is real [eq. (7)]. We also take into account the 
normalization condition of the states. In the second step we consider that all eigenstates 
contribute similarly and replace each | p by its average value l/Ne- We end up with an 
simple expression in the qubit (independent of the initial condition) and the form factor 
[eq. (7)] of the environment. We are also assuming that the level separation in //q is A so 
Eo = -El = A/2. 

Studying the average over initial conditions in the environment leads to the observa- 
tion that (Ai) = (A2) = 0(1/A'e). The remaining term can be computed along the same 
lines {i.e. same techniques and assumptions), however its calculation is both cumber- 
some and boring. We leave it to the enthusiastic student. Its result is 

Re(A3)=Re/^2(T-TO{l-gi(0i,02) + [l-^2(0i,02)]cos[A(T-T')]} (51) 

where gi {61,62) and g2( 01,^2) are geometric factors that depend on the initial state. 
Their values are given by 

8i{9u92)=g{ei)[\-g{e2)]+g{e2)[\-g{9i)] (52) 

g2{eu 02) = 2[1 -g{ei)] -g(02)[l -2g(0i)], (53) 

where g{0) = cos^ 9 + sin"* 0. 

The general solution for purity using this parametrization is 

P{t) = l-2?i^ f dx f dx'[\ + 8{x'/xn)-b2{x'/xn)] 
Jo Jo 

X [gi(0i,02)+g2(0i,e2)cosAT']+O(A4). (54) 

Two meaningful limits can be studied, the degenerate limit, in which the spacing of 
the internal Hamiltonian is much smaller than the inverse of the Heisenberg time of the 




FIGURE 8. (Color online) Numerical simulations for the average purity as a function of time in units 
of the Heisenberg time of the environment (spectator configuration, GUE case). For the coupling strength 
X ~ 0.03 we show the dependence of (P(f )) on the level splitting A in Hi and on the initial degree of 
entanglement between the two qubits (in all cases = 7i/A): 0i = (separable states), A = 8 (black 
triangles); 0i = ;r/4 (Bell states), A = 8 (blue rhombus); 0i = 0, A = (purple squares); 9\ ~ ;r/4, A = 
(red circles). We also show the corresponding linear response results (dashed lines) and exponentiated 
linear response results (solid lines). They are plotted with the same color, as the respective numerical data. 
In all cases A^e= 1024. 



environment (A ^ 1 / Th), and the fast limit (A 1 / Th) in which the opposite condition 
is met. In the degenerate limit purity decay is given by 

PD(O = l-A2(2-g0jr(f), (55) 

where 

2 , 

r(?) =?max{?,TH} + T — mm{t,Xn} . (56) 
3th 

The result is independent of 02 since a degenerate Hamiltonian is, in this context, 
equivalent to no Hamiltonian at all. The 0i -dependence in this formula shows that an 
entangled qubit pair is more susceptible to decoherence than a separable one. In the fast 
limit we obtain 

Pp{t) = 1 - [gi(0i, 02)KO +2THg2(01, 02)?] . (57) 

For initial states chosen as eigenstates of //q we find linear decay of purity both below 
and above Heisenberg time. In order for Pq to be an eigenstate of //q it must, first of all, 
be a pure state (in .^). Therefore this behavior can only occur if 0i = or 9\ = itjl. 
Apart from that particular case, we observe in both limits, the fast as well as the 
degenerate limit, the characteristic linear/quadratic behavior before/after the Heisenberg 
time similar to the behavior of fidelity decay. 

In Fig. 8 we show numerical simulations for the average purity decay. We average 
over 30 different Hamiltonians each probed with 45 different initial conditions. We 



contrast Bell states (62 = 7t/4, Oi = 7t/4) with separable states (62 = 7r/4, 61 = 0), 
and also systems with a large level splitting (A = 8) in the first qubit with systems 
having a degenerate Hamiltonian (A = 0). The results presented in this figure show 
that entanglement generally enhances decoherence. This can be anticipated since for 
fixed O2, increasing the value of 61 (and hence entanglement) increases both 62) 
and §2(^1 5 ^2)- At the same time we find again that increasing A tends to reduce the 
rate of decoherence, while a strict inequality only holds among the two limiting cases 
(just as in the one qubit case). From ^2(61, ^2) = 2 — gi{0i,02) — g{9i), it follows that 
{Pp — Pq)/X^ = g2{0i, O2) [r{t) — 2tTii] > 0. Therefore, for fixed initial conditions and t 
greater than 0, Ppit) > Poit)- This is the second aspect illustrated in fig. 8. 

The linear response formulae derived throughout this section are expected to work 
for high purities or, equivalently, when the second order Born expansion of the echo 
operator is a good approximation; it is of obvious interest to extend them to cover a 
larger range. The extension of fidelity linear response formulae has, in some cases, been 
done with some effort using super- symmetry techniques. This has been possible partly 
due to the simple structure of the fidelity amplitude, but trying to use this approach for 
a more complicated object such as purity seems to be out of reach for the time being. 
Exponentiating the formulae obtained from the linear response formalism has proven 
to be in good agreement with the exact super- symmetric and/or numerical results for 
fidelity, if the perturbation is not too big. The exponentiation of the linear response 
formulae for purity can be compared with Monte Carlo simulations in order to prove its 
validity. We now explain the details required to implement this procedure. 

Given a linear response formula /^r(?) for which P{0) = 1 [like eqs. (54), (55), or 
(57)], and an expected asymptotic value for infinite time Poo, the exponentiation reads 



This particular form guaranties that Pelr equals /^r for short times, and that 
limf_>oo/tLR(0 = The particular value of Poo will depend on the physical situa- 
tion; in our case it will depend on the configuration and on the initial conditions. 
Consider the spectator configuration. Let us write the initial condition as 



for some rotated basis. We assumed that for long enough times, the totally depolar- 
izing channel will act (recall that the totally depolarizing channel is defined as 
(fd[p] = ll/tr 11 for any density matrix p) on the coupled qubit. Hence, for the spectator 
configuration, the final value of purity is assumed to be 



PELR{t)=Poo + {l-Poo)exp 



i-/1r(0 



(58) 



l-P 



V/c(0)) =cos0i|OqOs)-fsin0i|lqls) 



(59) 




(60) 



Good agreement is found with Monte Carlo simulations for moderate and strong cou- 
plings, see fig 8. For weak coupling Poo depends on the coupling strength. 
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FIGURE 9. (Color online) We present the concurrence-purity CP plane. The area allowed for physical 
states is the gray area plus the set {{Q,P),P G [1/4,1/3]}. The Werner curve eq. (62) is shown as a 
thick black solid line. We show curves {{C{t)) ,{P{t))) as obtained from numerical simulations of the 
spectator Hamiltonian. We choose A = 1 and vary dim(J^) for two different coupling strengths X = 0.02 
(thick lines) and A = 0. 14 (thin lines). The resulting curves are plotted with different colors, according to 
dim{Jfg) as indicated in the figure. 



4.3. Two qubit concurrence decay 

In the previous section we studied purity decay of n qubits. Purity measures entan- 
glement with the environment, but we can wonder how decoherence affects the inter- 
nal quantum mechanical properties of the central system. Possibly the most important 
quantum mechanical property of a multi-particle system is its internal entanglement. We 
consider the simples scenario in which it is possible to have entanglement: a two qubit 
subsystem. Two qubit entanglement can measured with the concurrence C, a measure 
used extensively in the literature [49] that is straightforward to calculate. It is closely 
related to the entanglement of formation, which measures the minimum number of Bell 
pairs needed to create an ensemble of pure states representing the state to be studied [50]. 
Given a density matrix p representing the state of two qubits, concurrence is defined as 

C(p) =max{0,Ai-A2 -A3- A4} (61) 

where A,- are the eigenvalues of the matrix p {Oy Oy) p* {Oy (E) Oy) in non-increasing 
order. The superscript * denotes complex conjugation in the computational basis and Oy 
is a Pauli matrix [51]. The concurrence has a maximum value of one for Bell states, and 
a minimum value of zero for separable states. Furthermore, it is invariant under bi-local 
unitary operations. However, since concurrence is defined in terms of the eigenvalues of 
a hermitian 4 x 4-matrix, an analytical treatment, even in linear response approximation, 
is much more involved than in the case of purity, and has not been performed. 

We shall first explore a relation (first found in [52], partly explained in [53], and 
further studied in [17, 18]) between concurrence and purity. We show that this relation 
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FIGURE 10. (Color online) We show D [eq. (63)] which measures a "distance" in the CP plane between 
numerical curves for the RMT model and the Werner curve. In (a) its dependence on the size of the 
environment is studied for two cases. For A = 0.02 a finite value is approached, whereas for X =0.14 
algebraic decay is seen. The thick line, which is proportional to N^^, is meant as a guide to the eye. In (b) 
we plot D as a function of the coupling A, for various values of N^. In the large environment limit, and for 
/I <C 0.1, we have a noticeable deviation from the Werner curve. In all cases A = 1. 



is valid in a wide parameter range. Combining it with the appropriate formula for 
purity decay, we obtain an analytic prediction for concurrence decay. We compare 
our prediction with Monte Carlo simulations. It is essential that the two qubits do not 
interact; otherwise the coupling between the qubits would act as an additional sink (or 
source) for internal entanglement - a complication we wish to avoid. 

We study the relation between concurrence and purity using the CP-plane, where 
we plot concurrence against purity with time as a parameter. This plane is plotted in 
fig. 9. The gray area indicates the region of physically admissible states [54]. The upper 
boundary of this region is given by the maximally entangled mixed states. The Werner 

states pw = + ( 1 — a) |Bell) (Bell| , < a < 1 , define a smooth curve on the CP-plane 
(black solid line). The analytic form of this curve is [17] 



and will be referred to as the Werner curve. However, note that states mapped to the 
Werner curve are not necessarily Werner states. We perform numerical simulations in 
the spectator configuration. We compute the average concurrence for a given interval of 
purity using 15 different Hamiltonians and 15 different initial states for each Hamilto- 
nian. We fix the level splitting in the coupled qubit to A = 1 and consider two different 
values A = 0.02 and 0. 14 for the coupling to the environment. Fig. 9 shows the resulting 
CP-curves for different dimensions of Observe that for both values of A, the curves 
converge to a certain limiting curve as dim(^) tends to infinity. While for A =0.02, 
this curve is at a finite distance of the Werner curve, for A = 0. 14 it practically coincides 
with Cw(P)- Varying the configuration, the coupling strength, the level splitting, or the 
ensemble (GOE/GUE), gives the same qualitative results in the CP plane, for large di- 
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FIGURE 11. We show numerical simulations for the average concurrence as a function of time when 
both qubits are coupled with equal strength X to the environment. We consider small couplings A = 0.01 
which lead to the Gaussian regime for purity decay. The large dots show the result without internal 
dynamics, whereas the small ones are obtained for fast internal Rabi oscillations (A = 10 for both qubits). 
The theoretical expectation for concurrence decay based on eq. (66) are also plotted. The inset displays 
the corresponding evolution in the CP-plane with the same symbols as used in the main graph. In all cases 
A^e = 1024 




mensions. In some cases we have an accumulation towards the Werner curve, in others 
there is a small variation. 

To study this situation in more detail, consider a CP-curve Cnuml^)^ obtained from 
our numerical simulations, and define its "distance" D to the Werner curve as 



D= r dP|Cnum(i')-Cw(P)|. 



(63) 



The behavior of D as a function of the size of the system is shown in Fig. 10(a). For 
X =0.14 (dark dots), the error goes to zero in an algebraic fashion, at least in the range 
studied. In fact, from a comparison with the black solid line we may conclude that the 
deviation D is inversely proportional to the dimension of J^^. By contrast, for X — 0.02 
(light red dots), D tends to a finite value, in line with the assertion that the numerical 
results converge to a different curve. In Fig. 10(b) we plot the error D as a function of 
X, for different dimensions of J^. The results suggest an exponential decay of D with 
the coupling strength. The simplest dependence of the deviation in agreement with these 
two observations is 

We also plot the curves corresponding to this ansatz in fig. 10. Good agreement is 
observed for D <^ 1 . Notice the exponential decay of D with respect to X . One can thus, 
in an excellent approximation for large A, say that for large dimensions the limiting 
curve is the Werner curve. For A = 0, all studied couplings numerical convergence to the 
Werner curve was observed in the large A^e limit. 

Sufficiently close to P = 1, the above arguments imply a one to one correspondence 
between purity and concurrence, which simply reads C ~ P. This allows to write an 



approximate expression for the behavior of concurrence as a function of time 

Clr{t)=PLR{t), (65) 

using the appropriate linear response result for purity decay. The corresponding expres- 
sions for purity decay have been discussed in detail above. Eq. (65) has similar limits of 
validity as the linear response result for the purity. As it implies this approximation, we 
call it the linear response expression for concurrence decay. 

In those cases where the deviation from the Werner curve [see eq. (62)] is small and 
where the exponentiated linear response expression holds for the average purity, we can 
write down a phenomenological formula for concurrence decay, which is valid over the 
whole range of the decay 

Celr(0=Cw(WO)- (66) 

In Fig. 1 1 we show random matrix simulations when both qubits are coupled with equal 
strength A to an environment. We consider small couplings Ai = A2 = 0.01 which lead 
to the Gaussian regime for purity decay. We find good agreement with the prediction of 
eq. (66), except when we switch-on a fast internal dynamics in both qubits (A = 10) and 
consequently D is large. See the inset in fig. 1 1 . 

A similar analysis can be carried if one starts with a non-Bell state. However, as the 
reader can suspect this introduces more difficulties and more complicated expressions. 
The details can be found in [18]. 

5. CONCLUSIONS AND OUTLOOK 

The stability of quantum computation is the core of this course. We have given models 
for both fidelity decay and decoherence, based on random matrix theory. The most 
common exponential decay is a standard limit in both cases, but we have shown, that 
RMT not only reproduces the usual results in these cases, but provides a richer spectrum 
of behavior. Basically it introduces the Heisenberg time (of the total system or of 
the environment) as an additional time scale. This time scale in the case of fidelity 
decay defines the transition between quadratic (Gaussian) or perturbative and linear 
(exponential) or Fermi golden rule decay. Both regimes were known but with RMT 
the transition is well described. Furthermore it marks a clear difference to the behavior 
of integrable systems, where the size of the effective Hilbert space tends to be much 
smaller. This leads to faster decay, as the onset of quadratic behavior is much earlier. 

This brings us to an interesting discussion of great practical importance. Chaotic 
perturbations seem to be less harmful than integrable ones. This statement has given 
rise to many discussions with experimentalists [55], who, with good reason, say that 
the argument is doubtful because, with some adequate pulse in quantum optics or 
NMR (nuclear magnetic resonance), an integrable perturbation can be corrected for. 
Discussions of such corrections are found, e.g., in [56, 57]. 

Yet there is an important alternate aspect: By chaotizing an algorithm (after doing all 
the corrections possible) one can reach additional stabilization. This was first discussed 
for a quantum Fourier transform by Prosen and Znidaric [31]. They showed, that by 
introducing additional gates they can decrease correlations between perturbations and by 



consequence slow down fidelity decay. The idea to use random matrix theory to optimize 
such a method [58] took hold later and opens up very attractive alternatives. The basis of 
the idea is to randomly change the computational basis between gates. This field is just 
developing, and we refer the reader to [59] for recent advances. Fidelity freeze provides 
an additional option to improve fidelity decay. If we can show that pulse manipulation 
can eliminate, or reduce, the diagonal part of the perturbation, a freeze situation can be 
induced. The fact that even in the absence of strong spectral correlations a freeze can be 
typical, if not the ensemble average, opens a new alley of research. 
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